Long Tailed Maps as a Representation of Mixed Mode Oscillatory Systems 
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Mixed mode oscillatory (MMO) systems are known to exhibit some generic features such as 
the reversal of period doubling sequences and crossover to period adding sequences as bifurcation 
parameters are varied. In addition, they exhibit a nearly one dimensional unimodal Poincare map 
with a longtail. We recover these common features from a general class of two parameter family 
of one dimensional maps with a unique critical point that satisfy a few general constraints that 
determine the nature of the map. We derive scaling laws that determine the parameter widths of 
the dominant windows of periodic orbits sandwiched between two successive states of RL'' sequence. 
An example of a two parameter map with a unique critical point is introduced to verify the analytical 
results. 

PACS numbers: 82.40Bj, 05.45Ac 



I. INTRODUCTION 

Dynamical systems with disparate time scales for 
the participating modes often exhibit periodic states 
characterised by a combination of relatively large ampli- 
tude and nearly harmonic small amplitude oscillations. 
Such periodic states, called the mixed mode oscilla- 
tions (MMOs), and the associated complex bifurcation 
sequences have been observed in models and experi- 
ments in the area of chemical kinetics Q, |^ H, ^ |^, 
electrochemical reactions 0, S Q) biological systems 
, and in many physical systems 0, . These 
MMO systems typically exhibit the following features - 
(a) period doubling (PD) sequences and their reversal 
in multi-parameter space with respect to a primary 
bifurcation parameter keeping other parameters fixed, 
and (b) crossover to bifurcation sequences of alternate 
periodic-chaotic windows when one of the secondary 
bifurcation parameters is varied wherein dominant 
windows of periodicity increase in an arithmetical order 
which we refer to as period adding (PA) sequences 
iiiiillllll. 



periodic orbits of 1" kind constitute the period adding 
sequence. (Here, we use the conventional notation of 
to represent period orbits in MMO systems where L and 
s correspond to large and small amplitude oscillations 
respectively.) In the parameter space, wherein the 
reversal of period doubling sequences occur, periodic 
adding sequences are finite. Periodic orbits of 1^ type 
loose their stability in a period doubling bifurcation as 
the (primary) control parameter is increased, restabilize 
through a reverse period doubling bifurcation, and are 
eventually annihilated in a fold bifurcation [12| . In this 
parameter range, we have also studied the next maximal 
amplitude (NMA) maps obtained by plotting one 
maxima of the evolution of the fast variable with the 
next HI. ( The NMA maps can be regarded as a specific 
form of the Poincare maps.) These NMA maps show 
a near one dimensional unimodal nature with features 
of sharp maximum and a long tail. The sharpness of 
the NMA maps increases with the primary bifurcation 
parameter for a fixed value of the secondary bifurcation 
parameter. (See Fig. and 13). Also, see [H 
Several other MMO systems also exhibit similar features. 



As an illustration of the generic features of MMO 
systems, we collect some relevant results from our 
earlier study [13 . IT^ on Ananthakrishna's model (AK 
model) for a type of plastic instability j^. As in 
other MMO systems, this model also involves disparate 
time scales. The bifurcation portraits ( with respect 
to a primary bifurcation parameter, e) of the model 
show period doubling sequences and their reversal, 
gradually changing over to period adding sequences as 
the secondary bifurcation parameter, m, is increased. ( 
See Fig. QJ^jb.) As can be seen in Fig. the dominant 



A well studied example of MMO systems is the 
Belusov Zhabotansky (BZ) reaction system. Exhaustive 
theoretical/experimental studies for the BZ systems in 
two parameter space have shown that the NMA maps 
of these systems have a unimodal structure with a long 
tail and show a similar trend as a function of the control 
parameters as in the case of AK model. Other 

MMO systems which display similar features include 
lasers with a saturable absorber, autocatalytic systems 
and a number of oscillating chemical reaction systems 
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FIG. 1: Bifurcation diagrams for the AK model for a plas- 
tic instability with primary bifurcation parameter, e and sec- 
ondary bifurcation parameter, m. (a) Cascading period dou- 
bling bifurcations and their reversals forming a bubble struc- 
ture for m = 2.16, and (b) period adding sequences for 
m = 1.2. Chaotic regions exists in a vanishingly small pa- 
rameter regions sandwiched between successive l** periodic 
states. 



The periodic states representative of MMO systems 
have been shown to be associated vi^ith systems ex- 
hibiting global bifurcations in the form of approach 
to homocUnicity [l^, [H IHIH. Although, a variety 
of tools are available for the study of bifurcations in 
dynamical systems, understanding these features of 



MMOs has not been easy mainly due to the lack of 
adequate analytical tools to handle the global nature 
of bifurcations underlying the MMOs. However, since 
information regarding the nature and properties of 
periodic orbits of continuous flow systems are known to 
be embedded in Poincare maps, they have been used 
effectively to understand the underlying bifurcation 
mechanisms jl.6, 20, 24, 25] . For instance, considerable 
insight has been obtained through the Poincare maps 
derived from approach to homoclinicity [l^, |^ 112, ■ 
Since, most of the MMO systems exhibit high dissi- 
pation, there have been attempts to understand the 
dynamical behaviour of MMO system (specifically BZ 
system) by modelling its NMA map s as one dimensional 
discrete dynamical systems [HIM I27ll2ll29l|. 

The two distinct bifurcation features, namely, the 
reversal of bifurcation sequences and the period adding 
sequences have been modelled separately using one di- 
mensional maps. For instance, of the many mechanisms 
suggested for the reversal of period doubling sequences, 
a simple one involves maps having a single critical point 
along with nonmonotonous dependence of the control 
parameter jSOj . Another one is a unimodal map having 
a negative Schwarzian derivative 31j . Yet another one 
involves a competition between more than one critical 
point in the dynamics of the map [s^, 113, 13 ■ has also 
been shown that the transmutation of the U-sequences 
to the Farey sequences or specifically, its subset, the 
period adding sequences can arise in maps with more 
than one critical point with multiple parameters |35j. 
However, to the best of author's knowledge, none of 
these mechanisms apply to experimental and model 
MMO systems where one finds Poincare maps with 
unique critical point as well as smooth crossover from 
period doubling and its reversal to Farey sequences for 
a monotonous change in a secondary bifurcation param- 
eter. There has been no explanation for this crossover 
as well. Here, we show that a general class of maps with 
a unique critical point satisfying some broad constraints 
can show a gradual change from bifurcation diagrams 
having the reversal of the period doubling sequences to pe- 
riod adding sequences for a smooth change of parameters. 

To achieve this objective, we first abstract common 
features of MMO systems, and then impose them in the 
form of some general constraints on one dimensional 
maps with two bifurcation parameters in an effort to 
understand the above dynamical features of the MMO 
systems. ( Note that almost all MMO systems have two 
or more parameters.) Using these general constraints, 
we first analyse the basic mechanism responsible for the 
reversal of periodic sequences of RL'^ type which corre- 
spond to the dominant V sequence in MMO systems. 
(Admittedly, this is a much less difficult issue than the 
origin of reversal of all the bifurcation sequences in these 
maps.) 
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FIG. 2: Poincare maps from the AK model for m — 1.8 and 
(a) e = 190.6, and (b) e = 202.7. (Multiple folds arise as a 
result of finite dissipation.) 



We also show that the windows of symbohc sequence 
RL^ are the dominant windows in the parameter space 
and derive scahng relations for the onset of windows of 
periodic orbits sandwiched between successive windows 
of RL'^ sequence. The so derived scaling relations relate 
the slopes near the critical point and unstable fixed point, 
and are unlike the usual scaling relations derived in terms 
of bifurcation parameters. Investigation of the reversal of 
PD sequences in turn allows us to understand the mech- 
anism for the onset of the period adding sequences as 
well. Our analysis shows that the eigen value of the un- 
stable fixed point of the map plays a fundamental role in 



FIG. 3: Typical structure of the map-L described in appendix 
for /X = 4.0 and ^ = 18.0. The period one fixed point(p), 
critical point c and the first iterate /(c) are also indicated. 



controlling the window widths of the periodic orbits as a 
function of the control parameter in the family of long 
tailed maps. Hence, the geometrical shape of the map 
determined by the nature of unstable fixed point and the 
long tail structure is shown to be the underlying factor 
for the period adding bifurcation sequences. (For this 
reason, we refer to these general class of maps as long 
tailed maps.) 

We verify these scaling relations by constructing an 
example of a two parameter family of one dimensional 
maps with a geometrical shape similar to the NMA maps 
of MMO systems having a long tail. This example map, 
referred to as map-L, exhibits a smooth transition from 
period doubling bifurcation sequences and their reversals 
to period adding sequences for a variation of a secondary 
bifurcation parameter (Appendix). We emphasise that 
the example map-L is primarily introduced as an aid to 
verify the analytical results. However, to pin down the 
ideas, it will be useful to refer to the geometrical shape 
of map-L (Fig. |2J) which is seen to be similar to Fig. ^ 
Finally, we discuss the results and the correspondence 
with the continuous flow dynamical systems having a 
finite dissipation. 



II. DYNAMICS OF THE LONG TAILED MAPS 

In this section, , we first attempt to extract conditions 
that needs to imposed on one dimensional maps by 
studying MMO systems with generic features mentioned 
in Introduction section without any reference to any 
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functional form for the map. These conditions, stated 
below, are motivated by the study of the long tailed 
maps of MMO systems and are further aided by our 
detailed studies on the AK model 0, . Before doing 
so, we first attempt to capture the connection between 
periodic orbits of MMO systems and one dimensional 
maps. For the sake concreteness, we use Shilnikov 
scenario where the homoclinic contact is with respect to 
a saddle focus and motivate the issues using known facts 
of the AK model studied in detail [Hill Hi. 

We first identify the equilibrium saddle fixed point of 
a MMO system, as in the case of the AK model, with the 
origin of the map. In many of these experimental and 
model systems, interesting dynamics arises in the region 
between the first Hopf bifurcation resulting from desta- 
bilization of the fixed point and the restabilization of the 
fixed point from the periodic orbit through a reverse Hopf 
bifurcation. This is identified with back-to-back Hopf bi- 
furcation as in the case of the AK model also [^iTL IllllSq . 
The principal periodic orbit (PPO) born out of a Hopf 
bifurcation is identified fixed point p of the map ( see for 
instance Fig. 3). Further, careful observation of the bi- 
furcation diagrams of MMO systems show that the com- 
plex bifurcation sequence occurs in the region of the pa- 
rameter where the amplitude of the PPO decreases and 
approaches the reverse Hopf bifurcation point as the bi- 
furcation parameter is increased 0, 0, 

hUhJ. This im- 
plies that the fixed point moves toward the origin of the 
map. As the complex dynamics is terminated by a re- 
verse Hopf bifurcation as the parameter is increased, the 
height of the map also decreases as a function of a param- 
eter which we identify as the primary control parameter, 
u. These two statements can be written as two conditions 



fdf , ^df 



(1) 

(2) 



where c is the maximal point of the map. We expect that 
these conditions are satisfied by either of the two param- 
eters of the two parameter family of maps. Without loss 
of generality, we identify the parameter, fi, in Eq. ^and 
121 as the primary bifurcation parameter. (We shall use ^ 
for the secondary bifurcation parameter.) Here, we note 
that Eq. [21 implies that the map becomes increasingly 
sharp as fi is increased. As a consequence, the map also 
develops a long tail. Since, in most cases, the variables of 
continuous time systems are positive, we have taken the 
NMA maps to be limited to positive values. In Shilnikov 
scenario, as the orbit is reinjected close to the origin, the 
orbit spiral out along the unstable manifold. Closer the 
reinjection to the saddle focus, larger are the number of 
small amplitude oscillations. This can be captured by en- 
suring an intermittency channel in the one dimensional 
map. To mimic reinjection of the phase space orbit close 
to the saddle focus fixed point, we require = a 



and f{x,iJ,,^) a' < a for large x, where a and a' are 
non zero positive constants. Further, since we are inter- 
ested in maps with unique critical point c, there are two 
monotonic branches to the left and right of c. Further, 
as the whole map is coming down as a function of /i, the 
existence of a intermittency channel to left of c, should 
be anticipated. (Such a channel is actually observed in 
many systems lilllili^a^-) This channel is ensured 
if we impose 

/'(0,M,C)<1 (3) 
/"(0,Ai,C)>0. (4) 

The above conditions (Eqs. I1I4|I are equivalent to 
the observed features (a) to (c) of MMO systems as 
mentioned in the appendix. We shall, henceforth, refer 
to the general class of maps determined by these few 
constraints ( Eqs. liyj) as long tailed maps. 

To fix up the preliminaries, we start with some 
notations and definition. For the sake of brevity, 
wherever necessary, we will use f{x) in place 
f{x;fi,^). The n— th iterate of x is denoted by 
/"(a;) = /(/""^(a;)) = / o f'-^ix). The iterates of 
f{x;^i,^) are given by 9 : {a;o, a^i, a;2, • • • , x„, • • • } 
with the respective neighbourhoods given by 
{/o, /i, /2, • • • , In, ■ ■■}, where Xn 6 /«• A point q 
is fc-periodic if f''{q;f^,£,) = q, while /*(<;; 1 for 
0<z<fc — 1. In particular, we denote the period one 
fixed point by p. The eigen value of a fc-period cycle 

q is ^7(9) — \. If |A| < 1, then q is stable. Further, 
we will always choose the initial point, xq, to be the 
iterate that is closest to the critical point in the stable 
periodic window. Finally, we will deal with the non- 
trivial dynamics of the map in the interval / = [ci , co] , 
where /(c) = co and /^(c) — ci. Hereafter, we refer 
to the parameter windows of stable periodic orbits 
using the corresponding symbolic dynamics notation. 
(In this notation, to conform with the conventional 
representation, the symbol corresponding to the initial 
iterate xo has been dropped. For example, we represent 
the entire period three window by simply RL.) 

A rough understanding of the reversal of the period 
doubling sequences can be obtained by examining the 
changes in the stability of p. For a unimodular map, 
Eq.® implies that as fi is increased, as the fixed point 
p moves towards the origin, it loses stability leading to 
the first period doubling (PD) bifurcation, say at fj, — 
fi2- On the other hand, Eq.(^ implies the map as a 
whole is coming down which in turn means that p moves 
towards the critical point, c, as a function of /z. Once p 
crosses the point of inflection, ^{p) increases eventually 
culminating in a stable monoperiodic orbit when ^ (p) > 
— 1, say at /i = ^2 . Thus, if there are period doubling 
sequences beyond fi — 112, then their reversal has to occur 
within the parameter range /12 and /Xg. 

Since the map is comes down as a function of fi, there 
are two possibilities for reinjection, namely /^(c) falls to 
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Using the chain rule for iterates of maps, we have 





FIG. 4: Bifurcation diagrams for the Map-L. (a) Cascading 
period doubling with a bubble structure for ^ — 11.415, and 
(b) Bubble structure with dominant period three window for 
^ = 11.43. 



right of a point x = Xt at which f (xt, fi, £,) = 1 and 
/^(c) falls to its left. The latter clearly corresponds to 
the realization of period adding sequences. Any further 
increase in the parameter creates two more fixed points 
on the ascending branch of the map. 



A. Reversal of bifurcation sequences 

Before we proceed further, we give below some deriva- 
tives which will be used to establish subsequent results. 



df, 



dx 



(5) 



Using Eq. ij^l, we can write 



dx"^ 



'or 



dx 



n T^iPMm 



3=0 



Similarly, the parameter dependence of the iterates lead 
to 



ix).{7) 



In order to understand the mechanism of reversal of pe- 
riod doubling sequences exhibited by the unimodal long 
tailed maps, we restrict our discussion to a simpler prob- 
lem namely, the reversal of period doubling of the domi- 
nant sequences of the type RL'^ . (We will soon argue that 
these type of periodic orbits are the dominant ones.) 

To pin down the ideas, consider the map-L which 
satisfies the requirement of the long tailed maps. The 
dominance of the RL'^ periodic orbits are well visualised 
in the bifurcation diagrams of map-L (Fig. A 
prominent feature of the bifurcation diagram is the pres- 
ence of relatively dark bands of iterates ( corresponding 
to superstable orbits) that run across the bifurcation 
diagram connecting successive RL'^ periodic windows. 
(Similar features are routinely observed in bifurcation 
diagrams of MMO systems also.) These bands arise as 
a consequence of the increased stability of the iterates 
closer to the critical point of the map. It is clear that 
these bands share a similar sequential structure of 
iterates as the RL'' sequence, specifically the superstable 
point of the periodic orbit lies on the envelope of these 
bands. In the following analysis and later, we utilise this 
contiguous nature of envelope of bands to determine the 
parametric dependence of iterates. 

Consider the n iterates of A; + 2 periodic cycle which 
consists of one visit to the right of the critical point, c, 
and k visits to the monotonically increasing arm of the 
map. We examine the reversal of the sign of Eq.Q as 
fi is increased for periodic orbits of the form RL'^ as the 
reversal of the periodic sequence is assured if -g—{xo) 



changes sign at /i 
second iterate ( n = 



Mr 



To do this, first consider the 



2) of Eq.lTJ. Using xq — c, we get 



(8) 



Equation Q implies that the first term on the RHS 
of Eq.® is negative while the second term is positive 
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FIG. 5: The bifurcation diagrams for Map-L. (a) period 
adding sequence coexisting with reversal of period doubhng 
sequences at ^ = 12.8, and (b) period adding sequence with 
fold bifurcation at the accumulation point of RL" periodic 
windows for ^ — 18.0. Note the similarity with Fig. 0b. 



since f^(/(c)) and f^(/(c)) are negative. However, 

initially, lf^(/(c))| is small as J{c) refers to the 
tail part, and hence there is a region of /j. for which 
If > (c)| and thus, LHS of Eq.® is 

negative. As /i is increased ( as /(c) moves up towards 
c), |f|(/(c))| increases and ^^(c) turns positive, say 
at fi — fir2- (Note that the value of fXr2 is in principle 
a function of ^. For example, in the case of map-L, 
for ^ = 11.0, ^-(c) changes sign at ~ 2.2.) Beyond 



/i = fir2i-^{c) increases with /i. 

Now, consider evaluating Eq.(jZj) for successive n > 2. 
Consider, 

For periodic orbits of type RL'^, (/^C*^)) ^^'^ from 
Eq. O, f (/^(c)) < 0. Further, from the previous dis- 
cussion, since (c) is zero at ^ = fir^ : LHS of Eq. Q 

negative at this value. As /i is increased, ^^(c) turns 
positive and the second term on the RHS of Eq.@ is 
an increasing function of fi. Thus, the LHS which is 
negative for fj, < fj,^^ becomes positive at some value 

fi — firs > ^ir2 ■ In addition, note that | (c) | > | (c) | 
around ^ = /i^., . Similar arguments can be used to show 
that ^^(c), which is negative for /i < ^r„_.i will change 
sign at ^ — fir„ ■ Thus, the reversal of successive periodic 
orbits occur at increasing values of parameter given by 

(Note that these results are largely a consequence of 
Eqs.(l) and (2).) These results can be verified from 
the bifurcation portraits of map-L shown in Fig. 5a. 
Changes in sign of ^^(c), manifest as changes in 
the sign of the slope of the envelope of superstable 
iterates with parameter. It is clear that lower outermost 
envelope corresponding to the third iterate of critical 
point changes the sign of slope first, signalling the onset 
of reversal of the periodic sequence. 

A natural consequence of the reversal of periodic orbits 
is the presence of isolated bifurcation curves or the isolas. 
In the following, we trace the origin of of isolas using the 
first two constraints (Eqs. I1I2|) on the general class of 
maps. Again, we restrict our attention to periodic orbits 
of the form RL^. Consider using xo = c in Eq.® for an 
allowed fc-periodic RL^~^ , we get 

n—1 

n— 2 

We first note that excepting f^(/(c)) and |^(c), all 
other terms arc positive. Thus, as the visits are to 
the left of c, the sign of (c) is positive for all 
allowed values of fi. This positive curvature of the 
map and its parametric dependence, ^^(c) < for 
/i < firk, implies that a fc-periodic orbit is created due 
to a tangent bifurcation as fj, is increased. For further 
increase in ^, f^{c) comes down, triggering a period 
doubling bifurcation making the periodic orbit unstable. 
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Increasing beyond /i^^. reverses the period doubling 
since ^^(c) becomes positive, ie., the region around 

c of f^{x) moves up eventuahy restabihzing the fc— th 
periodic orbit. Further increase in ^ destroys the stable 
periodic orbit in another tangent bifurcation (see Fig. 
5b). The creation of a periodic orbit, its destabiliza- 
tion and restabilization, followed by its destruction 
constitutes an isolated bifurcation curve (isola). This 
happens for each of the allowed periodic orbit of the 
form RL^ . From these arguments, it is clear that the 
isola corresponding to the smallest allowed period orbit 
contains all isolas of the allowed periodic orbits of larger 
k, with the smallest one being the maximum periodic 
orbit. In addition, there exists a maximum periodic 
orbit of RL^ that is allowed. In the case of map-L, it 
is clear that the value at which different iterates show 
reversal in sign ^ (c) are different as can be seen from 
Fig. 5a and the reversal of all the allowed periodic 
orbits of RL'^ type occurs beyond ^ ^ 4.05 (^ — 8.0) 
and the maximum allowed periodic sequence period five 
{ k — 3). Indeed, similar isola features are well doc- 
umented in bifurcation diagrams of MMO systems H [13. 

It is possible to extend these arguments to some 
other sequences of periodic orbits as well. For example, 
consider periodic orbits of RL^R^"^ type. Similar 

arguments as before show that — '9^2 (c) is positive 



and 



(c) changes sign from negative to positive, 
hence forming isola structures. However, there are other 
kinds of periodic orbits in the allowed symbolic sequence 
(for example, RL'^R^™^^), whose reversal is interrupted 
by the homoclinic points. Extending the above type of 
arguments to such cases will depend on the nature of 
the sequences and appears difficult. 

As discussed in the introduction, one important fea- 
ture of MMO systems is the dominance of the parameter 
windows of stable periodic orbits of RL'' type. Using 
Ea. (|10|l . we show that the width of the parameter win- 
dows of these stable periodic orbits depend critically on 
the structure of the map as determined by Eqs.QJ and 
In particular, we will show that the ratio of the 
widths of successive periodic windows of the RL'^ type is 
decreasing. Consider Eq. (|10|l for k and k + 1 stable peri- 
odic orbits. In comparison with the fc-th periodic orbit, 
there is an extra term in the product arising from an it- 
erate falling on the monotonically increasing arm. Since 
the fi value of the {k + l)-th periodic orbit is larger than 
that for fc-th one, which in turn implies /(c) has moved 
towards c, we have 



n9^(rw)<n9^(/"w)- (10) 

n=2 n=2 

Further, since /(c) moves towards the critical point, 

^(/(c))| increases. Noting that ^^(c) is negative, 
whose magnitude increases with fc ( or /x), it implies 



that f''~^^{x) exhibits a sharper minimum around c 
than that for the fc-th periodic orbit. Assuming that 
the extent of the neighbourhood around /'^(c) can be 



approximated by i^* 



c)^/2, the window 



width of the fc-th periodic orbit can be approximated by 
the product of F'= and ^(c). Since |F'=| < \F''+^\ and 



< 



-(c) 



it clear that the window width 



of the fc-th periodic orbit is larger than that of (fc + 1) 
periodic orbit. Since the value at which the reversal 
of PD sequence of the type RL'^ occurs depend on the 
secondary bifurcation parameter ^, it is possible that 

there may not be a change in the sign of ^^(c) beyond 

a certain value of fc. ( Note that if does not change 
sign, then, there cannot be a change in sign for to > fc. 
) Thus, only a few periodic sates of RL'^ would be seen. 
However, for the period adding sequences, the number of 
states of RL'' can go unbounded. For this case, a channel 
at small x should open up. The onset of a channel can 
be anticipated based on Eq. ^ which implies that the 
graph of the map comes down with fi, and /(O, ^, /i) — a, 
i.e., nonzero value as the origin of the map. A power 
series representation for the monotonically increasing 
part with conditions /'(0,/i, ^) < 1 and /"(O,^,^) > 
ensures that there is one point to left of c that has 
/'(a;,/i,^) = 1 for all values of the parameter. Then, for 
RL'^ periodic sequences, the window widths of periodic 
orbits of the type RL'^ decreases slowly. 



B. Scaling relations for periodic orbits 

Having argued for the predominance of the periodic 
orbits of sequence of RL'^ type in the parameter space, 
we show that periodic windows contained in the chaotic 
region bounded by the windows of periodic orbits de- 
noted by RL'' and RL^^^ have much smaller parameter 
width compared to that of RL^ and RL^^^ . Numerical 
results on the map - L are also presented in support of 
the analytical result. 

The PA sequence also manifest in NS-unimodal maps 
as a part of MSS sequence 38] . In the MSS sequence, 
between any successive RL^~^ and RL^ windows, there 
exist windows of allowed sequences of the type RL''[S], 
where [5'] is a sequence consistent with the allowed sym- 
bol sequences [39^. i.e.. 



RL 



<---'<RL^[S] <RL^. 



Considering only the increasing behaviour of long tailed 
maps, we assume that the MSS ordering corresponding to 
the periodic orbits of NS-unimodal maps may be applied 
for these types of maps also. 

Let the bifurcation values jik correspond to the last 
generation fc-period orbit in the MSS sequence, i.e., for 
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FIG. 6: (a) Map-L and the periodic orbit of form /?L^ and 
Rli^R for ^ = 3.7831 and p. = 3.99 respectively = 18.0). 
(b) Verification of the scaling of f^(a;o'^) with |Ap|~^ for 
RL^R periodic sequences obtained from the example Map- 
L. The dashed line is drawn for a linear fit. 



responding to the simplest symbolic sequence in the re- 
gions i?/ where stable periodic windows are of the type 
KL^R\S\. Within this sequence, the most dominant ( 
in the sense of the MSS sequence) periodic sequence is 
RL^R which has one period higher than the basic se- 
quence RL^ . Using Eq.©, for a /c— th periodic orbit ( 
/'^(xo) = xq), the eigen value of the periodic orbit, A, is 
given by. 



and for stability 
point, we have Ap 



n— fc— 1 

^ n 

n=0 

< 1.0. 



dl 

dx 



{Xn) 



(11) 



(For the period one fixed 
< 1.0 for ^ < /i2.) Within 



the first order, we approximate all iterates faUing in the 
neighbourhood J7 of p as having the same slope as at p. 
In the case of long tailed maps, owing to the sharp nature 



of longtailed maps, 



{P) 



is large, for ji > ji^. For the 

sake of clarity, consider the sequence RL^R which is part 
of the sequence RL^R\S] in i?/. Clearly, the structure of 
the map ensures that one iterate of this sequence belongs 
to U (see for example map-L in Fig. 6a ). In the case of 
RL^R sequence, f^^'^{xQ) falls in neighbourhood U of p. 
Using Eq. jSJ for stability of the periodic orbit. 



df, .df 



n— 1 



< 1.0 



(12) 



where Xq is an initial point chosen to be in the neighbour- 
hood of the critical point for reasons that will become 
clear soon. Since long tailed maps are sharply peaked ( 
Eqs.lQJ and |(5J ), almost all the iterates are trapped in 
the intermittency channel like region near the origin of 
the map. (This point is best illustrated by considering 
the concrete example of map-L shown in Fig. 6a, where 
RL^ and RL^R are shown.) Thus, for large fc, noting that 
iterates other than xq and Xk+2 change only marginally, 
for stability of the periodic orbit, we have 



C 



^(.o)^(x.+2) 



< 1.0 



/i > /ife no more cycles of period k are present. Thus, the 
parameter values are ordered in the following way: 

Ai3 = As < ^3 < /*4 < /i4 < ^5 < • • • 

where jl\ refers to degenerate homoclinic bifurcation 
po| (wherein f^{c) — p) of the periodic point p between 
k-th and (fc -|- l)-th periodic orbits. 

Using a degenerate homoclinic bifurcation as the 
boundary, the parameter space for all the periodic win- 
dows between RL^~^ and RL^ ( of period A; -1-1 and k + 2) 
can be divided into two regions, as Ri : [/ifc+i, and 
Rii : /ife+2]. Consider the periodic window cor- 



where C ( = Y\n=i f'{xn)) is a constant factor. Since 
Xk+2 falls in the neighbourhood of p S C/, 
|Ap|(> 1.0). Hence, 



Of 



iixk+2, 



dl 

dx 



(2^0) 



< [c\xA\ 



(13) 



Recall that for any periodic orbit, one iterate is always 
located in the neighbourhood of the critical point in the 



interval 



^PD ^SN 



'-0 







where 



and Xq correspond to 



the onset of the periodic orbit through a saddle node bi- 
furcation at ^ and the destabilization of the periodic 
orbit through a period doubling bifurcation at fJ,!^^ ■ As 
the parameter is increased from /i = , the iterate xq 
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can be regarded as traversing from Xq^ to . Then, 
Eq. (14) implies that for the stabihty of the periodic 
orbit RL'^R, the iterate (.tq) is restricted to a smaller ex- 
tent of the neighbourhood of the critical point than that 
corresponding to RL'^ in order to compensate for the in- 



crease of 



11(^/^+2) 



|Ap|. Assuming a uniform change 



in xq with respect to the parameter /i as it traverses from 



to Xq^ and using the fact that Ap 



■\p 3> 1.0, the pa- 
rameter windows of RL'^R is smaller by a factor Ap than 
that for RL'^ . Similar arguments can be used to show 
that any sequence with 2m + 1 extra visits to the right 
than RL'' would have a window width smaller by a factor 
^-(2m+i)^ (The assumption here is that these visits to 
the right are close to p.) Since any RL'^RIS] sequence 
consists of 2m -|- 1 right visits over and above the ba- 
sic sequence RL'^ , the width of this periodic window is 
vanishingly small for large to. Thus, 



a/ 

dx 



< c- 



,_(2m+l) 



(14) 



More specifically, Eqs. (13) and (14) allows us to 
get a scaling relation between the onset of the periodic 
windows of type which can be verified by 

using the map-L. For this we fix the control parameter, fi 
at the fold bifurcation points of the periodic windows so 
that the equality holds in Eq. (14). Under this condition. 



Ap — C 



Eq.(14) ( for TO = 0) implies that f^C^^o ) 
which is almost independent of k in RL^R. To verify 



this, we have plotted 



which exhibits linear scaling behaviour for the onset 
values of RL"R windows using the map-L (Fig. 6b). In 
a similar fashion, we have shown a plot of ^{xq^) as a 
function of jApj-^^'^+i) for various m values correspond- 
ing to periodic orbits of RLF'R^"^^^ type (Fig. 7). Note 
that, even though large k approximation is used, the 
scaling relations work well when k is not large (Figs. 6b 
and 7). 

In the region Rn , the allowed periodic orbits are made 
up of the sequence RL'^R^IS] which has at least two more 
iterates to the right than the corresponding RL'' orbit. 
(See Fig. 8a.) The smallest allowed periodic orbit in this 
region is RL'^R? which has exactly two iterates more than 
the periodic orbit RL'^ both of which lie in the neighbour- 
hood U ofp. ( Note that even though the actual distance 
of these iterates may be substantially away from p, due 
to sharpness of the map, the slopes at these points can be 
approximated by the slope at p.) Thus, approximating 
the derivatives of the map for these two iterates falling 
in the neighbourhood U of p to be Ap, the eigen value of 
the periodic orbit RL^R^ is seen to increase by a factor 
Ap compared to RL'^ orbit. Hence, we have 



as a function of lAr, 



|Ap| 



(15) 




FIG. 7: Verification of tlie scaling of ^{x^^) with 
|^^|_(2^+i) periodic sequences obtained from 

example Map-L. 



Following arguments presented for Rj, we see that the 
iterate Xq is restricted to a smaller extent around the 
critical point by a factor A^ ^ which in turn leads to the 
width of RL'^R^ being smaller by the same factor com- 
pared to RL'^. Using the equality, the above relation has 
been verified numerically for the map-L with the control 
parameter, /x, kept at the fold bifurcation point. Figure 
8b shows the scaling for the onset of the periodic or- 
bit RL'^R'^ with Ap ^ which confirms the analytical result. 

It is to be noted that in the chaotic region between RL'^ 
and RL'^^^, there exist stable windows of other symbolic 
sequences also. In the above, the scaling is shown only 
for a particular sequence of periodic windows, namely 
RL'^R"^. However, the structure of the map ensures that 
any periodic orbit with periodicity higher than that of 
RL'^R and RL^R^ in the regions i?/ and Rn respec- 
tively, typically have more number of iterates falling into 
the neighbourhood U of p. This structural feature along 
with the large value of Ap ensures that the other window 
widths of these periodic orbits are increasingly small. / / 



III. CONCLUDING COMMENTS 

One of main motivating factor for the study of a 
specific type of Poincare maps, namely the NMA maps, 
of MMO systems is that these maps constitute the 
dimensionally reduced form of the continuous time 
system; thus these maps are expected to retain all the 
key features of the MMO systems. It is well known 
that a dominant feature of homoclinic bifurcation to an 
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FIG. 8: (a) Periodic orbit sequence of RL*R and RL-^R^ for 
11 = 3.990345 and fi = 4.001927 for ^ = 18.0 for the Map-L 
(the maps are indistinguishable from each other on this scale), 
(b) Verifications of the scaling relation ^{xq'^) with |Ap|~'^ 
for RL"R? periodic sequences obtained from Map-L. 



unstable limit cycle (Gavrilov-Shilnikov scenario) ^41] is 
the existence of isola structure. This is observed in the 
region of periodic-chaotic sequences in MMO systems, 
and even in those cases which show incomplete approach 
to homoclinicity 0, 0, ^^l- ( See below.) Since the 
isola structure is preserved at the level of long tailed 
NMA maps, as a part of the study of long tailed maps, 
we have undertaken to investigate the origin of these 
isolas by studying the reversal of bifurcation sequences 
(within the limited scope of RL^ type of orbits). For 
example, the RL*' sequences were shown to be created 
and subsequently destroyed through a fold bifurcation 
leading to isola structures. Hence, in the parameter 
region where period adding seguences ( of type RL*' ) 



manifest, these isolas are prominent and an isola with 
larger periodicity is embedded within the isola with 
a lower periodicity leading to a progressively ordered 
structure of isolas. This feature mimics the ordered 
structure of the isolas present in certain types of MMO 
systems. 

From the point of view of continuous time systems, 
the reversal of period doubling sequences is a direct 
consequence of the back-to-back Hopf bifurcations in 
these systems 0, ^3, ^2 ■ For a two parameter dynam- 
ical system, a degenerate Hopf bifurcation gives rise to 
back-to-back Hopf bifurcations and the existence of the 
periodic orbits is confined to this interval in the pa- 
rameter space. Within the region of back-to-back Hopf 
bifurcations, any bifurcation wherein periodic orbits are 
created has to be matched by a reverse bifurcation of the 
same kind wherein the periodic orbits will be destroyed. 
In other words, in these continuous flow systems also, 
the periodic orbits born in a fold bifurcation vanishes 
in another fold bifurcation creating an isola structure. 
Moreover, the amplitude of the period one orbit arising 
from the first Hopf bifurcation approaching the other 
Hopf bifurcation (while going from the periodic state to a 
fixed point) gets progressively smaller. Since the period 
one orbit is equivalent to a fixed point of the NMA map, 
this feature translates to the condition ^(p) < that 
has been imposed for the changes in the structure of the 
map. A consequence of this condition is that the first 
period doubling occurs due to the changes in the shape of 
the map in the neighbourhood of period one fixed point, 
while the rest of the period doubling bifurcations are 
due to the presence of the critical point. This difference 
in the cause of the first period doubling bifurcation 
and higher order bifurcations is refiected in the first 
Feigenbaum exponent bein g di stinctly different from the 
rest of the PD bifurcations Eg . Equivalently, a large 
period two window is seen as a distinct feature of the al- 
ternate periodic-chaotic sequences of MMO systems also. 

We have shown that the increasing sharpness of 
the map and the consequent long tail nature, as the 
parameter is tuned, determines the dominance of the 
RL'^ periodic states. Heuristically, the dominance of 
the RL'' type periodic orbits can be related to the 
arguments on dissipation in the following manner. The 
primary features of the structure of the map relevant 
to the predominance of period adding sequences is the 
large eigen value (in absolute terms) for the period one 
fixed point, the long tail structure, and the presence 
of an intermittency channel. In these class of maps, 
apart from the neighbourhood of the critical point, there 
exists a large tail region wherein the slope is small. By 
definition, these regions are related to high dissipation, 
since the spread of the iterates in its next iteration is 
suppressed. (Given an neighbourhood Ii in the tail 
region, I2 = f{Ii;ii,^) I\. ) One type of sequence 
of periodic orbits which include visits to the regions of 
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high dissipation with no visits to the neighbourhood 
containing p are the RL'' sequences. High dissipation 
favours periodic orbits, hence the stability of these 
periodic sequences are enhanced at the expense of the 
chaotic regions. In the same spirit, the windows of 
periodic orbits occurring within the chaotic region 
between any two successive RL'' sequences, involving at 
least one iterate in the neighbourhood of the unstable 
fixed point and favouring negative dissipation (local 
expansion) have window widths smaller by a factor of 
Xp. These analytical results for the onset of periodic 
orbits of the type RL'^R™' show that the slope at the 
point of onset of the periodic orbit, namely ^^(a^o) is 
related to the slope of the unstable fixed point Xp, and 
scales linearly with A~™, where m is the number of 
visits to the neighbourhood of the unstable period one 
fixed point. The relationship between the slopes near the 
critical point and the fixed point essentially controls the 
width of such periodic orbits. R must be emphasised that 
these scaling relations derived are very different from 
those usually derived where the widths of the periodic 
orbits are evaluated as a function of the parameters. 
Again, map-L has been used only to verify these scaling 
relations. The scaling relation suggests that the width 
of the parameter window is nearly independent of the 
visits to the left of critical point of the map and only the 
visits to the neighbourhood of p determine the width of 
these type {RL^) of periodic windows. For a map with 
a sharp maximum and a long tail, the eigen value of 
the unstable fixed point will be large and consequently 
smaller windows widths for the periodic orbits contained 
in the chaotic region sandwiched between RL^~^ and 
RL''. 

A necessary condition for the occurrence of the 
reversal of period doubling sequences in unimodal maps 
have been stated in the form of the first two conditions 
defined by Eqs. (1) and (2), which are constructed 
keeping in mind the changes in the structure of the 
NMA maps of the MMO systems. The correctness 
of these conditions is validated by the fact that the 
dominant bifurcation sequences of the long tailed maps 
and the MMO systems have similar features. Further, 
the similarity of the bifurcation diagrams obtained 
from the example map-L ( belonging to the class of 
long tailed maps ) is similar to the MMO systems and 
in particular bifurcation sequences from AK model. 
Since the dissipation involved in the higher dimensional 
continuous time is related to the separation of time 
scales operating in the system, a larger separation of 
time scales implies near one dimensional nature as well 
as the long tailed structure for the map. Since the 
scaling relations incorporate the essential ingredients of 
the structure of the map, namely, the property of long 
tail and the strength of the unstable fixed point {Xp), the 
scaling relations derived for the one-D maps can be used 
as an effective check for the kind of dynamics involved in 
the higher dimensional systems. These observations are 



subject to usual criticism applicable to maps considered 
as reduced dynamical systems. In higher dimensional 
systems, it is conceivable that the interaction between 
several parameters can lead to reversal of period dou- 
bling and period adding sequences. The present work 
only attempts project the entire complexity into a few 
constraints. Clearly, these constraints on the maps can 
be satisfied in variety of combinations of parameters in 
higher dimensional space. 

In summary, we have traced the dominance of period 
adding sequences as a part of complex alternate periodic- 
chaotic sequences in MMO systems to the long tailed na- 
ture of the associated Poincare maps. In the process, 
we have partially elucidated the origin of the reversal of 
the bifurcation sequences with specific reference to RL^ 
type of orbits which are the dominant sequences in MMO 
systems. The results are derived based on some general 
conditions alone without taking recourse to any particu- 
lar analytical form of the map; the map-L is used only 
to numerically confirm the analytical results. Further, 
the example map-L which satisfies the general constraints 
mimics the distinctive features of higher dimensional con- 
tinuous flow MMO systems, namely, the reversal of PD 
bifurcation sequences and the alternate periodic-chaotic 
sequences. The specific feature of the dominance of win- 
dows periodic orbits over the chaotic regions in the long 
tailed maps also captures the equivalent feature for the 
continuous time MMO systems. Since, many maps of ex- 
perimental MMO systems exhibit generic features similar 
to those used here, the scaling relations derived for the 
onset of the periodic orbits of the form RL^R™ can be 
taken to serve as a check for the correctness of the anal- 
ysis. In particular, in BZ systems, while earlier experi- 
ments do show the dominance of the period adding se- 
quences |29j , the existing results are not accurate enough 
to verify these scaling relations. However, we believe that 
careful experiments should bear out our results. Indeed, 
we expect that careful experiments on a number of other 
experimental systems wherein the geometrical shape of 
the Poincare map is similar to Fig. 2 [ESHEli could 
be used to validate our results. Another example would 
be the laser systems where accuracy of control on pa- 
rameters is generally considerably high. These results 
also indicate the ubiquity in the qualitative dynamical 
features of physical systems from widely differing origin, 
exhibiting alternate periodic-chaotic sequences. 



IV. APPENDIX 

In this appendix, we provide an example two param- 
eter map which has the following broad typical features 
observed by the NMA maps of MMO systems: a) non 
zero origin, b) a sharp symmetric maximum positioned 
asymmetrically, and c) asymptotic long tail to the right 
of the maximum. It is a generalisation of the Bountis 
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map|30j given by 

Xn+l = f{Xn;fJ.,0 

- /x) + r{x - cY 



f{x;fi,0 = a + 



with m > I and a > 0. We refer to this map given by 
Eq. (Al) as map-L. This map is used mainly to verify 
the analytical results derived in text. We shall use 
only two parameters /i and ^ of the several parameters 
l,m,a, r, si, S2 and c. The other parameters of the map 
/(a;;/i,^) are tuned such that the model map looks 
similar (in structure) to one dimensional maps of the 
MMO systems (Fig. 2). For our numerical study, 
we use the following values for the map parameters : 
I = m = 4.0, a = 3.0, r = 0.6, Si = 1.5, S2 = 0.12 and 
c = 7.0. We shall use /i and ^ as primary and secondary 
bifurcation parameters respectively. 

The bifurcation diagrams are constructed using /i as 
a control parameter for fixed value of ^. Fig. 3 shows a 
typical shape of the map defined by Eq. (Al). This map 
has a unique smooth maximum with inflection points 
on either side of the maximum. For properly chosen 
initial conditions, the asymptotic tail ensures that the 
nontrivial dynamics of the map is restricted to the 
interval [/^(c; ^, i^), /(c; /z, ^)]. It can be easily verified 
that the map has a negative Schwarzian derivative 
everywhere on the positive real line. Since a > 0, for an 
appropriate choice of a, there are no fixed points in the 
interval [0, c] for small values of ^. For a range of values 
of the parameter ^ < the map has only one fixed 
point, p, which is in the interval [c, /(c)]. For ^ > 
two more fixed points, p2 and pa, located in the interval 
[0, c], are born in a fold bifurcation at /i = /i*(f). ( Both 
p and P3 are unstable, and p2 is stable.) 

In the following, we briefiy describe the bifurcation se- 
quences of the map /(x; ^, with respect to the primary 
bifurcation parameter /i for various values of ^. For small 
but fixed value of ^, an increase in the value of /i leads 
to a decrease in the value of p. For ^ < 8.0, the unique 
fixed point p is stable for all values of /i. Beyond ^ = 8.0, 



the fixed point p loses stability in a period doubling bi- 
furcation as /i is increased which is regained in a period 
undoubling bifurcation thus forming a bubble structure 
(Fig. 4a). For further increase in ^, the number of nested 
bubbles grc^s^is 2" with n = 1, 2, 3 • • • . For ^ - 11.9, 
a period three window opens up in the bifurcation dia- 
gram separating the period doubling sequences from the 
reverse period doubling sequences (Fig. 4b). Further in- 
crease in the value of ^ unveils period adding sequences 
of arithmetically increasing periodic windows coexisting 
with the bubble structure (Fig. 5a). This sequence of 
large parameter windows of stable periodic orbits (in /i) 
is found to be of RL'^ type sequences in the two letter 
symbolic dynamics language. For ^ — 12.8, we find that 
the reversal of the dominant bifurcation sequences occur 
beyond ^ ^ 4.05. At ^ = 18.0 (^* 13.0), as we increase 
fi, at fi ^ fi* = 4.8, the region of the map for small values 
of X makes contact with the bisector giving rise to a fold 
bifurcation with the creation of a pair of stable and un- 
stable fixed points, p2 and ps (Fig. 5b). Keeping ^ > 
as the parameter fj, is increased towards /it*, the chan- 
nel between the bisector and the map gets progressively 
narrow, thereby stabilising periodic orbits of higher peri- 
odicity (See Fig. 3 and also Figs. 6a, 8a). The positions 
of the iterates of higher periodic cycles change marginally 
to accommodate the next higher periodic cycle as addi- 
tional iterates are squeezed in this intermittent channel. 
Thus, the intermittency or the fold bifurcation point (fi*) 
is an accumulation point for the arithmetically increas- 
ing period adding sequences. Beyond the accumulation 
point, no further bifurcations are observed for the stable 
periodic point (pa). This map reproduces most bifurca- 
tions features of the AK model. ( C omp are bifurcation 
diagrams here with those in T^] and '48'].) 
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